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We study the motion of ballistic electrons in a single miniband of a semiconductor superlat- 
tice driven by a terahertz laser polarized along the growth direction. We work in the semiclassical 
balance-equation model, including different elastic and inelastic scattering rates, and incorporat- 
ing the self-consistent electric field generated by electron motion. We explore regions of complex 
■ dynamics, which can include chaotic behavior and symmetry-breaking. Finally, we estimate the 

' magnitude of the DC current and voltage that spontaneously appear in regions of broken-symmetry 

for parameters characteristic of real semiconductor superlattices. 
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I. INTRODUCTION 



q_) ' Since the seminal work of Esaki and Tsu [1] predicting negative differential conductivity in the current-voltage 
characteristics of a semiconductor superlattice (SSL), SSLs have been the subject of intense theoretical and experi- 
mental study. SSLs consist of thin layers of alternating semiconductors - typically 10-100 layers ranging from 10 to 
500 A in thickness - whose different band structures create a periodic potential along the growth axis, essentially 
forming a one-dimensional, artificial crystal. The superlattice period can be controlled during the growth process and 
is generally much larger than the period of the constituent materials. Consequently the conduction and valence bands 
break into narrow allowed and forbidden regions termed minibands, with typical widths ranging from 1 to 100 meV. 
For an introduction to SSL design and properties, see reference [2]. 

Recent experiments on SSLs in the presence of high-frequency helds have revealed many novel, nonlinear effects, 
including photon-assisted tunneling [3], a reduction of the current during laser irradiation [4], absolute negative 
conductivity [5] , resonant increases in the current [6] , and, most recently, deterministic chaos in the sequential resonant 
tunneling regime [7]. Our previous work [8] suggests that chaos should also occur in the miniband regime for a SSL 
driven by an AC field. In this brief note, we recall some of the principal results of that work and report our current 
studies of symmetry-breaking effects in electron transport. Using numerical simulations we establish the possibility of 
spontaneous DC current and voltage generation in response to a purely AC drive. Our results cover the whole range 
of frequency/amplitude of the applied field, extending a recent investigation of the high frequency, strong field limit 
[9]. 



II. BALANCE EQUATIONS 

We consider electrons subject to an external field E ext — Eq cos(fit) using the tight-binding dispersion relation, 
e(k) = y[l — cos(fca)], for SSL period a, minibandwidth A, electron energy e and crystal momentum k. We assume 
Q \ that the frequency of the external field is not large enough to support Zener tunneling to higher minibands, and that 
O ■ the electric field is uniform across the SSL. To model the electron dynamics, we generalize the semiclassical balance 
equations of Ignatov et al. [10] by including the self-consistent electric field resulting from electron motion, obtaining 
[8], 



V = -eE tot {t)/m(e) - lv V (1) 
e = -eE tot (t)V - j s (e - e ) (2) 
E sc = -47TJ*/ k ~ aE sc (3) 

where E to t = E ext + E sc is the total electric field, including the self-consistent contribution; j = —eNV is the current 
density; V(t) and e(t) are the average electron velocity and energy, respectively; 7„, 7 e , and a are phenomenological 
relaxation rates for the average velocity, average energy and self-consistent field, respectively; k is the average dielectric 
constant for the SSL; e n is the temperature-dependent equilibrium energy; N is the carrier concentration, assumed 
uniform across the SSL; and m(e) — pr^rs ' m ° = 'KS 1 > ls en?ec tive mass of the electrons. 
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The balance equations remain invariant under the symmetry transformation: 

t _> t + ^±It, y -> -v, £ sc -» -s sc , (4) 

where T is the period of the applied field and n — 0, 1,2, in other words, if {V(t), e(t), E sc (t)} satisfy (0)-(||), 
then so do {-V(t + ^-T), e{t + ^±lT), -E sc (t + ^±lT)}. Physically, this represents the lack of preferred spatial 
direction for a SSL without DC bias; under this transformation the electron velocity and the self-consistent field 
reverse direction, but the electron energy remains unaltered. 

Symmetric solutions also remain invariant under transformation (^); hence V(t + 2 "+ 1 = — V(t) and E sc (t + 
j~ T) = —E sc (t). For such cases the time-averages of both the velocity and the self-consistent field are zero. However 
there may also exist solutions which do not satisfy this symmetry, for which V and E sc have non-zero average values; 
these symmetry-broken solutions correspond to the spontaneous generation of a DC current and bias. 

In the limiting case j v — j e — 0, (|l|)-@ can be transformed to the equation describing the Resistively Shunted 
Junction (RSJ) model of a Josephson junction (JJ) [8,11] This model possesses a thoroughly studied symmetry 
analagous to (Q) [12,11]. When biased by an AC current of frequency f2, the JJ may develop a DC component to its 
voltage difference satisfying the phase-locking condition (U) oc nfi, where n is an integer, U is the voltage across the 
J J, and angular brackets indicate time-averages [12]. 

III. NUMERICAL STUDIES 

In our earlier work [8], we showed that the balance equations give rise to chaotic behavior for a wide range of 
damping and driving parameters. Readers should consult this reference for a full description of the numerics and 
scalings used in the simulations. Briefly, however, there are three important frequencies: the drive frequency, tt; the 
Stark frequency of the driving field, u> s — eaEo/h; and the frequency of the cooperative oscillations associated with 

the self-consistent field, lue — \J We 2 N a 2 A/ nh 2 which is usually in the terahertz range. As it depends only on SSL 
material parameters, we use loe as the scale of inverse time, and throughout the rest of this article we consider the 
dimensionless frequencies D,/loe, ^s/^e, "Jv/^e, Je/^e and cx/uje without explicitly showing the scaling. 

As in many studies of the J J [12,11], we examine the dynamics for fixed damping while varying the parameters 
of the applied field. Figures la and lb indicate the region of chaos and symmetry-breaking, respectively. Note the 
similar general shape of these regions; further work shows that for greater damping the region of chaos diminishes 
while the region of symmetry-breaking retains this general shape [8,13]. 

Fig. lc plots the time-averaged self-consistent field, {eaE sc /hujE}, as a function of external frequency, uj s , for the 
frequency Q = 0.9. Phase-locked symmetry-breaking where (eaE sc /?uvE) ~ nfl/ujE for n = ±1, ±2 is clearly visible. 
A careful examination reveals that (eaE sc /?uvE) also depends weakly on lo s ; therefore the phase-locking relation is 
only approximately satisfied. 

From (^J) it immediately follows that (V) oc (E sc ) for non-chaotic solutions. Consequently symmetry-breaking should 
be experimentally observable as a DC current and bias across an unbiased SSL illuminated by a terahertz laser. For 
representative SSL parameters of a = lOOl, ]V = 3x 10 15 cm" 3 , A = 22me^, and e = 13 [6], u E = 3.3 x 10 12 s _1 . 
Taking (eaE sc /TnoE)max = 1-5 yields a DC bias per SSL period of 3mV. Furthermore, a lateral area of 8/im 2 would 
result in a current of lfiA. Both these magnitudes should be experimentally detectable. 

IV. CONCLUSION 

We have studied the complex dynamics of semiclassical balance equations that model the motion of electrons in the 
miniband of a semiconductor superlattice exposed to an external ac field. The equations include the self-consistent 
field created by electron motion. Numerical simulations indicate the possibility of broken-symmetry solutions, which 
occur in roughly the same laser frequency and electric field ranges as chaotic dynamics. The symmetry-breaking 
engenders a DC bias and current that could be experimentally detectable. 

We would like to contrast this work with recent theoretical [14] and experimental [7] work predicting and observing 
chaos for SSLs in the incoherent tunneling regime when the electric field breaks into domains of different magnitudes 
[2]. Our model considers the coherent motion of electrons within a miniband in a spatially homogenous electric field. 
A complete theoretical understanding of when the applied field breaks into domains has not yet been attained, so it is 
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encouraging that chaos is predicted in quite different modeis applicable in different regimes of SSL design and applied 
field characteristics. 
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FIG. 1. For a weakly damped SSL, chaotic (dark symbols) and regular (white areas) regions (a); symmetry-broken (dark 
symbols) and symmetric (white areas) regions (b); time-averaged self-consistent field as a function of applied field strength (c) 
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